Thermal emissivity for finite three-dimensional photonic band gap crystals 
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We discuss the results of computer model for the thermal emissivity of a three-dimensional pho- 
tonic band gap (PBG) crystal, specifically an inverted opal structure. The thermal emittance for a 
range of frequencies and angles is calculated. 
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I. INTRODUCTION 

It has been shown that there exist periodic dielectric 
structures such that for a range of frequencies, no elec- 
tromagnetic waves can propagate in any direction for all 
polarizations |0, ^, ^, |4[ g. The properties of such an 
infinite crystal have been worked out in detail ^, 
Such work is typically limited to non-absorptive mate- 
rials, since otherwise all waves become evanescent. 

Much theoretical work has been done on the proper- 
ties of finite one dimensional photonic band gap (PBG) 
crystals js), including recent calculations of the thermal 
emissivity of such one-dimensional structures . The 

strong angular dependence of the gap effect with a one- 
dimensional structure has motivated successful experi- 
mental work with three-dimensional structures |pd]| . The 
three-dimensional model we will use here is an inverted 
opal structure, and with the use of a computer program 
published by Stefanou et al. O, 13 , the transmittance, 



reflectance, and absorbance of a finite number of layers 
of the crystal is calculated, as well as the complex band 
structure of the corresponding infinite crystal. As we 
shall show, our structure exhibits a broad band gap in 
the emissivity that has virtually no angular dependence, 
and as such, also suppresses the thermal emittance of 
the crystal in an omnidirectional fashion. We also pre- 
dict emissivity enhancement at the photonic band edge. 

For a given crystal structure, we can calculate T and 
TZ, the transmission and reflection coefficients of the crys- 
tal as a function of the frequency m and the wave vector 
k of the incident wave. If we consider a crystal with real 
indices of refraction, then TZ + "T = 1 from conservation 
of energy. However, if there is an imaginary part to the 
any of the indices of refraction, there will be a non-zero 
absorbtion coefficient A — 1 — TZ — T . Calculation of 
this quantity is useful because of a Kirchoff 's second law 
which states that a material's thermal emmittance £ is 
proportional to its absorbance (and that they are equal 



for a blackbody). One can obtain the absorbance by 
simply multiplying Planck's blackbody power spectrum 
by the absorbtion coefficient. So by modelling the ab- 
sorbtion of several layers of a three-dimensional crystal 
coated on the surface of a thick, highly absorptive sub- 
strate, we are also modelling the thermal emissivity E of 
the entire structure. The three-dimensional superlayer 
acts as a passive filter, which alters the emissivity of the 
underlying graybody substrate. 



II. THEORY 
A. One-dimensional case 

We employ software developed by Stefanou et al. |l^, 
that uses a generalization of the matrix-transfer tech- 
niques developed for one-dimensional crystals, such as 
thin-film dielectric stacks. In the one-dimensional case, 
the T and TZ are and |rp respectively, where t and 
r are the complex transmittance and reflectance coeffi- 
cients, respectively. These coefficients are linearly related 
to one another by a 2 x 2 matrix M, where 



and 
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The column vectors on each side of Eq. (Q) repre- 
sent the fields outside each side of the crystal. Then 
t — 1/Afii and r — M21/M11. This matrix M can be con- 
structed for a one-dimensional structure simply by multi- 
plying the transfer matrices for the individual layers. For 
a homogenous plate, these can be readily computed by 
taking into account the the refiection and transmission 
at the interface at the edges of the plate with the Fresnel 
equations, along with the accumulated phase (and pos- 
sible change in amplitude if the index of refraction has 
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an imaginary component) from propagating though the 
material. 

Because the phase accumulated from the wave prop- 
agating though a given layer depends on the portion of 
the wave-vector that is perpendicular to the interfaces 
of the plates, the frequency lo at which the gap is cen- 
tered is highly dependent on the angle of incidence of 
the initial radiation. As a consequence, while the on-axis 
thermal emissivity at a given lo might be suppressed, off- 
axis emissivity might have a band-edge enhancement, as 
shown in earlier work |9[ p^ . 

B. Three-dimensional model 

The model used by the Stefanou and co-workers 
also allows for a finite number of parallel planes of 
spheres of differing index of refraction from the host ma- 
terial in which they are embedded. For a given frequency 
w, the incident plane wave is expanded into spherical 
waves about an origin centered at one sphere, and then 
the scattered wave from a single sphere in one plane 
is calculated by satisfying the boundary conditions for 
Maxwell's equations at the surface of a single sphere. The 
scattered wave for the entire infinite plane of spheres is 
just the sum of the scattered waves for a single sphere 
expanded about a translated origin, with a phase factor 
corresponding to the portion of the incident wave vector 
that is parallel to the plane of spheres. By summing this 
first-order scattering of all the spheres, except the sphere 
at the origin, and then expanding the result into spherical 
waves about the origin, we can get a better approxima- 
tion of the true wave that is incident on the sphere at the 
origin. 

Here, the restriction that each plane of spheres has the 
same two-dimensional periodicity allows for the difficult 
part of this expansion to be calculated only once per 
frequency, which is what makes this method much more 
numerically feasible than other approaches. 

By solving the boundary conditions again, with the 
spherical wave expansion of the incident plane wave plus 
this first order scattered wave, we solve for a second-order 
scattering off the sphere at the origin. By symmetry, all 
the spheres of the infinite plane have this same scattered 
wave. Hence, the total scattered wave for the plane is 
again the sum of the scattered waves for a single sphere, 
expanded about a translated origin, with a phase factor. 
The end result is then expanded into plane waves and 
added to the incident wave to to solve for the total trans- 
mitted and reflected field. By repeating this process for 
each plane of spheres, the total field for the entire crystal 
is obtained. 



III. RESULTS 

The three-dimensional crystal chosen for our calcula- 
tions is a diamond fee inverted opal (Fig. |l|) with the 




FIG. 1: A unit cell of the inverted opal structure used. 
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FIG. 2: The angular dependence of the emissivity for the 
inverted opal structure with four periods. The frequency is 
given in units of c/a, where c is the vacuum speed of light and 
a is the length of the unit cell in Fig. 

complex relative dielectric constant ei = 12 + i/10 for 
the space around the spheres and 62 = 1 inside the 
spheres, which are of radius 0.30618621 a. Four layers 
of the unit cell are stacked on a highly absorptive sub- 
strate (e = 12 + 7i) which is eight orders of magnitude 
thicker than the height of the total crystal, so that it can 
be considered essentially infintely thick. The emissivity 
dependence on normalized frequency and polar angle for 
both polarizations is given in Fig. ^ 

Several important features are noted. Firstly, there is 
a large band-gap in the emissivity that is largely angle in- 
dependent. Secondly, the emissivity of the graybody sub- 
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FIG. 3: Emissivity for a one-dimensional crystal with 16 pe- 
riods at several angles from normal. Contrasted with a three- 
dimensional crystal, it has a much stronger angular depen- 
dence. 
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FIG. 4: Part of the band structure (left) and transmission 
(right) of the infinite, non-absorbing crystal with ei — 22. 

strate goes exactly to unity at the band-edge resonance 
frequencies. This means that the emission is enhanced 
up to the pure blackbody limit rate at these frequen- 
cies. This observation is consistent with our previous 
one-dimensional calculations |l^, and recent experi- 
mental observations | pd| . As noted in our previous work, 
around the band-edge resonances the three-dimensional 
photonic crystal superstrate acts as an antireflective coat- 
ing, allowing all incident radiations at the frequency to 
be absorbed — which is the definition of a blackbody. 

You will notice that a complete gap does not exist for 
this structure, as is the case for most inverted opals, as 



you can see emission is not suppressed for a particular an- 
gle for the center gap frequency of 2.27 c/a. While other 
structures might show more promise of having complete 
stop bands, opals arc unique in that large crystals can be 
grown by self-assembly of individual spherical particles, 
which can then be doped out to provide the inverted opal 
structure we have presented. 

A complementary set of plots for a one-dimensional 
crystal is given in Fig. ^. For the one-dimensional crystal, 
indices of refraction and thicknesses chosen such that the 
absorption gap is centered at the same frequency and of 
the same width as the three-dimensional crystal. This is 
achieved by setting ei = 2.6, di = 0.6 a, 62 = 1.44, and 
^2 — 0.81 a where di and ^2 are the thicknesses of the the 
two layers and ei and 62 are the corresponding relative 
dielectric constants. Here, a is once again the length of 
the edge of the unit cell of the three-dimensional crystal, 
so all frequencies are comparable. 

While the frequencies corresponding to suppressed 
emission in the one-dimensional crystal are highly de- 
pendent on angle, the emissive gap effect in the three- 
dimensional crystal is nearly independent of angle. 
Though, as mentioned before, a band structure cannot be 
calculated for a finite, absorptive structure, one can cal- 
culate the band structure for the infinitely extended crys- 
tal with the absorption removed by setting Im(ei) — 
and removing the backplane. However, the existence of 
the gap effect in our configuration is highly dependent on 
the boundary condition of the thick absorptive backplane 
and the presence of absorbtion in the crystal itself, and 
no gap effect is present if they are removed. This geom- 
etry does exhibit gap behavior at a different frequency, 
though, if ei is changed appropriately (Fig. 

IV. CONCLUSION 

We have calculated the thermal emissivity of a three- 
dimensional PBG crystal using extensions of tested nu- 
merical techniques developed for ID crystals. Our 
matrix-transfer approach is conceptually much simpler 
than the traditional QED method of solving the non- 
trivial boundary conditions involved, and also not nearly 
as computationally intensive as the Finite Difference 
Time Domain (FDTD) simulations used by many groups 
for solving Maxwell's equations in PBG crystals. Our 
technique also takes into account the essentially finite 
nature of the structure, crucial for experimental work 
where often it is possible to fabricate only a few layers 
of crystal. Our results agree qualitatively with recent 
experiments by S. Lin [pl| . 
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